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Abstract. We show how to construct families of stationary hydrodynami- 
S-J cal configurations that reproduce the observed vertical gradient of the rotation 

velocity of the extra-planar gas in spiral galaxies. We then present a simple 
| model for the lagging halo of the spiral galaxy NGC 891, which is in agreement 

with the Hi observations. Our method is based on well known properties of 
' baroclinic solutions, and it is an elementary application of a much more general 

and flexible method. 
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• i-H , 1. Introduction 

ft 

Observations at different wavelenghts show that spiral galaxies are surrounded 
by a gaseous halo. This extra-planar gas is mu ltiphase: it i s detected in H I (e.g., 
ISwaters. Sancisi. Sz van der Hulstl 11997ft Ha (iRandl l2000ft . and X-ray observa- 



ISwaters. Sancisi. &: van der HulstJ 11997ft Ha (IRandl [2000 ) , and X-ray observa- 
tions (|Wang et al.ll200ll : IStrickland et al.ll2004 ) . In particular, high-sensivity Hi 
obser vations of edge-on galaxies like NGC 891 ( Swaters et al.lll997l : lFraternali et al.l 
2004j|) and UGC 7321 dMatthews fc Woodll2008l 'l reveal neutral gas emission up 



to large distances from the plane and show the presence of a vertical gradient 
in the gas rotation velocity (see Fig. In addition, the study of face-on galax- 
ies has revealed the presence of vertical motions of neutra l gas often associated 
with holes in the disk Hi distribution dPuche et al.l ll992. Boomsma et al., this 
conference) . 

Clearly, the two major issues of the origin and of the dynamical state of 
the extra-planar gas are strictly related. For example, the halo gas could be the 
result of cosmological accretion (i.e., infall of extragalactic gas; Binney, this con- 
ference), or could have been ejecte d from the plane through a galactic fountain 
mechanism (Sha piro Sz Fieldlll976ft . or it could have had a mixed origin. One 
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might expect different structural and kinematical configurations for these cases. 
Here we focus on the problem of the dynamical state of the extraplanar gas. 

There are two "extreme" kinds of explanations for the extra-planar gas kine- 
matics (low rotation and vertical motions): the ballistic models and the fluid 
stationary models. Ballistic models describe the gas as an inhomogeneous col- 
lection of clouds in ballistic motion, subjected only to the gravitational potential 
of the galaxy. A well-known example of ballistic model is the galactic fountain, 
which describes how ionized gas is ejected from the galactic plane due to stellar 
winds an d supernoya exp losions and then cools and falls back ballistically onto 
the disk (iBregmanl 11980). Ballistic models are effective at explaining vertical 
motions of the cold (Hi ) and warm (Ha) gas components. However, the obser- 
vations indicate a considerable morphological and kinematical regularity of the 
extra -planar gas (see, e.g., the total Hi map of NGC 891 in iFraternali et ail 
2004 and Fig. Q] here) , which may be difficult to understand in the context of 
purely ballistic models (see also ICollins. Benjamin. &: Randll2002l ) . 

The observed regularity may be explained more satisfactorily in the context 
of fluid stationary models. In these models, the gas is taken to be a rotating 
fluid in stationary equilibrium, without motions along the R and z directions; all 
the thermodynamical quantities are therefore time- independent, and the galaxy 
gravitational field is balanced by the pressure gradient and the centrifugal force. 
Until now, this approach has not been fully explored in all i ts possibilities, and 
only a few attempts have been made (e.g., see lBeniaminl2002T ). Here we present a 
short account of our approach and of i ts main advantages, while a full discussion 
will be given in a forthcoming paper ([Barnabe et al.ir2004h . 



2. Stationary models: the standard approach 



The simplest fluid models are constructed by solving the stationary equations 
of hydrodynamics (written in the standard cylindrical coordinates R, z, and <p), 
under the assumptions that vr = 0, v z = 0, and v v = SIR (in other words, the 
system is in a state of permanent rotation): 
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+ n 2 R, 



where p, P and O are the density, the pressure and the angular velocity of the 
gas, and 3>tot is the total gravitational potential of the system. Owing to the 
axial simmetry, all the physical variables are functions of R and z only. 

A commonly adopted ap proach for the solution of Eqs . I|T|) is that based on 
the Poincare-Wavre theorem ( Lebovitzlll967l : lTassou]|ll980f ). According to this 
theorem, the effective gravitational field at the r.h.s. of Eqs. (JJ) can be obtained 
from an effective potential $ e fj if and only if O = £l(R). In this case, the gas 
density, pressure, and temperature are all stratified on ^eff; the relation between 



P and p is necessarily barotropic [i.e., P = P(p)], and 



j(.R). The last 



property, i.e. the so-called "cylindrical rotation", is in clear disagreement with 
the observed vertical gradient of the extra-planar gas rotation velocity in spiral 
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galaxies. Thus, in standard applications, one fixes $tot an d the radial trend of 
0(i?) (thus fixing ^eff)- Then, one solves the equation VP = — pV^eff f° r the 
density, with assigned boundary conditions and with specified P = P(p); for 
= one recovers the standard (barotropic) hydrostatic equilibrium solutions. 
Clearly, rotational velocities changing with z are excluded here, and this would 
seem to argue against the fluid stationary approach. 

However, in the next Section we will show that it is possible to construct 
baroclinic equilibrium solutions with a negative velocity gradient along z: this 
approach is a simple applicatio n of a more general technique that will be pre- 
sented bv lBarnabe et al.l (|2004T ). Note that baroclinic solutions have been stud- 
ied before for problems ran ging from geophysics to the theory of sunspo ts to 
galactic dynamics (e.g., see lRosselandlll92^ : lTassoullll980l : IWaxmanlll978l . and 
references therein). In particular, isotropic, self-gravitating axysimmetric galaxy 
models can be interpreted as baroclinic flu id configurations, showi ng streaming 
velocities often decreasing with z (e.g., see lLanzoni fc Ciottill2003h . 



3. Stationary models: baroclinic solutions 

At variance with the standard approach, here we start by assuming a gas density 
distribution p(R, z) (vanishing at infinity). Thus, for a given & tot (R, z), the first 
of Eqs. ((J) can be integrated to obtain the gas pressure 

P(R,z) = j™p^dz>. (2) 

In general, P can not be expressed as a function of p only, and so our system is 
barocline. Its temperature distribution is obtained from the perfect gas equation 
of state, 

T=^EL. (3) 

kp 

Inserting Eq. © in the second of Eqs. (^Q) and integrating by parts one obtains 
a remarkable "commutator-like" expression for the rotational velocity 



K _ f° f°p d ®tot dp a$tot 
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Clearly, due to the baroclinic nature of the solution, the quantity v v will depend 
on z. Unfortunately, the construction of barocline solutions is not as easy as it 
appears: in fact, not all the density distributions produce physically acceptable 
solutions, that is solutions for which v z, > everywhere. However, using Eq. Q 
we can state a few sufficient conditions that can be used as guidelines to choose 
the density distribution p and to obtain physically acceptable solutions. From 
now on, for simplicity we assume that the gas is not self-gravitating. Two 
sufficient conditions are the following: 

(i) when the total potential is homeoidally stratified with axial ratio </<j>, i.e. 

z 2 

*tot = *tat(*) s £ 2 cxR 2 + — , 0<g*<l, (5) 

4 
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(as in Binney's [1991] logarithmic potential and Evans' [1994] spheroidal 
potentials), if one assume a gas density distribution of the form 

z 2 

p(R,z) = Pl (R)p 2 (m), m 2 ocR 2 + —, < q„ < 1, (6) 

q 2 

a sufficient condition to have v 2 ^ > everywhere is 
dpi(R) dp 2 {m) 

^ir- ' ^r^ ' and (7) 

In particular, the third condition requires the gas density distribution to 
be stratified on flatter homeoids than the equipotentials. 

(ii) When the total potential is generated by a razor thin uniform disk, (i.e. 
<3? tot = 27tG'Eqz), a necessary and sufficient condition to obtain v 2 > 
everywhere is dp(R, z)/dR > 0. 

In particular, point (ii) suggests that if the stellar disk contribution to the 
total potential is dominant in the central regions, a "sufficient condition" to have 
physically acceptable solutions is to take a centrally depressed gas distribution. 
The physical reason for this condition is very simple. In fact, suppose to fix a 
generic p(R, z) distribution in a uniform (vertical) gravitational field. The gas 
distribution is not stratified on <J> tot and so must be rotating. According to the 
second of Eq. (^Q) this means that, at any fixed z, the gas pressure must increase 
with R. This request can be stated in terms of the gas column density pdz' . 
In fact from Eq. (|3J) it results that in a vertical gravitational field the square 
of the gas rotational velocity is just proportional to the radial gradient of the 
column density, that must be positive. It is worth mentioning that this trend is 
consistent w ith the observed Hi surface density distribution for several spiral 
galaxies fsee lCavatte et al.1ll994h . 



4. A simple model for the NGC 891 lagging halo 

As a simple application of baroclinic solutions we have attempted to reproduce 
the observed negative gradient of the rotational velocity of the extra-planar gas 
in NGC 891. The model here should be considered exploratory, in the sense that 
we are not interested at reproducing perfectly the kinematical data, but only to 
provide a reasonable, physically acceptable model for the extraplanar gas, and 
to show that baroclinic solutions deserve deeper investigations in the context of 
the extra-planar gas kinematic. 

We consider a very simple galaxy model that consists of two components, 
namely a razor-thin exponential stellar disk 

GM d f°° Jp(fej)ej*l ~_ R __ 

^ MZ) = -^J (l + k 2 f/ 2 R =R~ d > Z =R- d > (8) 



and a logarithmic dark matter halo (Bi nnevlll98lh with constant asymptotic 
velocity vq 



^ R = ^ 



ShaioCR, z) = ^ In [Rl + R 2 + ^ ) , R h = ^. (9) 
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The galaxy bulge is not considered here. In Eqs. (jHJ)-® and are the 
disk mass and scale radius, and % are the halo potential scale radius and 
flatte ning parameter, and J„ is the n-th order Bessel function of the first kind 
(e.g.. lBinnev fc Tremainelll987h . 
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Figure 1. Rotation curve of NGC 891 for the gas in the plane of galaxy 
(filled circles) and for a strip from 30" to 60", corresponding to about 2.25 
kpc, above and below the plane (open circles). The solid curve represents the 
circular velocity in the equatorial plane of the model galaxy, while the dashed 
curve is the rotational velocity of the baroclinic solution 2.25 kpc above the 
galactic plane. 



Following the simple criteria illustrated in points (i) and (ii), we adopted 
the trial function 

^ Z ) = ^U + I C " W - = & + (10) 

Rq (1 + m z )Pl z g| 

where po, i?o and h g are the central density, the radial and vertical scales of the 
gas distribution, and a, j3 and 7 are dimensionless constants. For a > 0, Eq. ()10j) 
describes a centrally depressed (i.e. toroidal) gas density distribution, and the 
exponential is a scale-height modulating function. Note that the equations for 
this mod el must be solv e d num erically, and the description of the code will be 
given in Ba rnabe et alJ (J200J). After a few attempts we fixed R& = 3 kpc, 
M d = 8 x 1O 1O M , v = 220 km s _1 , R h = 5, q h = 0.71, R = 2 kpc, h g = 3/2, 
a = 1, f3 = 4, 7 = 1/2 and q g = 1/4. For this choice the model circular velocity 
reproduces the observed rotation curve of NGC 891 (see Fig. ^ solid line), and 
%> > everywhere. According to Eqs. ©-©, the quantities v v and T do not 
depend on the value of pq. 

The meridional section of the isorotation surfaces are shown in Fig. [5] as 
solid lines: v v decreases with z over the main body of the distribution, with the 
exception of the region R < 2R& and z ~ 0, where v v varies in a non-monothonic 
way with z. 

The dashed curve in Fig. ^ corresponds to a horizontal cut of Fig. [21 at 
z ~ 2.25 kpc, superimposed to the observed values of the Hi rotation velocity 
at the same z. The match between the observed points and the model curve is 
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Figure 2. Meridional section of isorotation {solid, in km s 1 ) and isothermal 
(dashed, in Kelvins) surfaces for the model in Eqs. (|%|l- H10JI . 



surprisingly good, in particular when considering that the exploratory nature of 
the presented model. 

The model isothermal surfaces are also shown in Fig. |2] as dashed lines. 
As expected, the gas distribution is quite hot, with temperatures in the range 
10 5 < T < 10 6 K. The hotter region is "bubble" near the symmetry axis above 
(and below) the equatorial plane, while the gas t emper ature decreases steadily 
approaching the galactic disk. In lBarnabe et alJ ( 2004f ) we present and discuss 
(soft) X-ray total luminosity and surface brightness maps of our models. In order 
to compare the model with the observations the value for the central gas density, 
po, must be known. In fact, one could proceed as follows: while the dimension- 
less part of Eq. (|1U|) can be fixed to reproduce the rotational and temperature 
fields, the observed total X-ray luminosity could be used to determine po- With 
this parameter fixed, one can then discuss other, astrophysically relevant model 
properties (see next Section). 



5. Discussion and conclusions 



We have presented fluid models for the extra-planar gas in spiral galaxies. The 
approach is based on the class of hydrodynamical equilibria known as baroclinic 
solutions. In particular, we have showed how very simple baroclinic configu- 
rations can be caracterized by a decrease of rotational velocity with increasing 
z similar to that observed in the extra-planar gas of spiral galaxies. We re- 
mark that baroclinic solutions are very flexible: for instance, if we require v v 
to reach the systemic velocity for large \z\ (i.e., the halo becomes hydrostatic), 
one could repeat the previous analysis with a factorized gas density distribution 
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as p(R, z) = f(R, z)p e ($tot), where / ~ 1 at large \z\ and p e is an hydro- 
static equilibrium solution in the galactic potential. Even more general cases 
can be obtained by "perturbing" cylindrical rotation solution, i.e., by assuming 
Pe = Pe(&efi)- For such distributions several interest ing results similar to those 
reported in points (i) and (ii) can be easily derived (iBarnabe et al-lfeOQ^ b 

We conclude by discussing the major limitations and open problems related 
to the use of a fluid approach: 

1) Homogeneity. Observationally the gas is multiphase: what is the rela- 
tionship between the (hot) single-phase gas described by fluid models and the 
Hi and Ha components of extra-planar gas? 

2) Absence of vertical motions. Observations show vertical motions of the 
order of 50 — 100 km/s, while in the approach described in Sect. 3 such motions 
are excluded: is this enough reason to abandon stationary fluid solutions? 

About point 1), a first important insight will be given by the calculation, 
for a chosen baroclinic distribution, of the local and global cooling times. In fact, 
while the global cooling time sets the rate at which heat must be furnished to 
the distribution in order to assure stationarity, the local cooling time indicates 
where clouds will condense out of the smooth density distribution. It would be 
interesting to follow the trajectories of such falling clouds under the drag force 
of the gas and with initial conditions given by the baroclinic solution. Thermal 
conduction could be another important ingredient in the evolution of the clouds. 
Of similar interest would be the study of the interaction of cosmologically ac- 
creted cold gas with a hot rotating halo. In particular, it would be interesting to 
measure evaporation times and to find out whether and on what time-scales the 
kinematical regularity of the fluid solution is transferred to the infalling cloud 
population. 

As to point 2), in addition to the comments above, we also note that baro- 
clinic solutions can be generalized to include meridional motions. This general- 
ization is not trivial but can be obtained, at least in the geostrophic l imit (i.e., 
for small Rossby n umbers), by standard expansion techniques (e.g., IWaxmanl 
119781 : lTassouHll980h . Numerical methods are required to investigate existence 
and properties of solutions with faster meridional motions. In any case, it can- 
not be excluded that also meridional motions may play a role in the kinematics 
of the hot gas. 

Thus, it is clear that the ballistic and fluid approaches both give important 
(but restricted) informations about the state of the extra-planar gas. A bet- 
ter understanding of the problem will be obtained by the investigation of the 
interaction of clouds with baroclinic fluid configurations using time-dependent 
numerical hydro dynamical simulations. Numerical simulations will also clarify 
the fate of the "torus-like" structure common in baroclinic distributions. 
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